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Abstract 

We predict cascade hypernuclear production spectra expected in the forthcoming J-PARC experiment. In the Green's 
function method of the distorted wave impulse wave approximation with the local optimal Fermi averaging t-matrix, 
we can describe the H~ production spectra in the continuum and bound state region reasonably well. Predictions to 
the high resonlution spectra at J-PARC suggest that we should observe E~ bound state peak structure in (K~, K + ) 
spectra in light nuclear targets such as 12 C and 27 Al. 
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1. Introduction 

Investigation of nuclear systems with strangeness 
opens up a way to understand dense matter such as 
the neutron star core. Since strange quarks are nega- 
tively charged and cancel the positive proton charge, 
they are favored in charge neutral dense matter. As 
a result, many models predict hyperons would ap- 
pear at around 2po ; and A may share a similar or 
larger fraction to neutrons at very high densities. 

In describing highly dense matter, we clearly need 
information on BB interaction not only for S = 
0,-1 but also for S < —2 such as EN, but spectro- 
scopic information on cascade (S) hypernuclear sys- 
tems are severely limited at present. While old emul- 
sion data suggest a deep S~ -nucleus potential (~ 

-24 McV) [1], a shallow potential ( 15 MeV) is 

suggested from the twin hypernuclear event found in 
a nuclear emulsion [2] . This shallow potential is also 
supported by the distorted wave impulse approxi- 
mation (DWIA) analysis of S~ production spectra 
in the bound state region [3,4]. 

In order to extract as much information as pos- 
sible from the data available at present, we need to 
investingate the E~ hypernuclear production spec- 



tra by (K~ , K + ) reaction on nuclear targets in both 
of the continuum as well as the bound state region. 
Since the DWIA analysis in [4] strongly rely on the 
absolute value of the inclusive E~ production yield, 
it is necessary to verify the consistency with the pro- 
duction spectra in the quasi free (QF) region [5]. For 
this purpose, the Green's function method of DWIA 
would be a useful tool, where the continuum and 
bound state spectra can be described on the same 
footing. 

In this paper, we investigate the S~-nucleus po- 
tential through S~ production spectra in the contin- 
uum and bound state region in the Green's function 
method of DWIA with the local optimal Fermi av- 
eraging t-matrix (LOFAt), in which the S~-nucleus 
potential effects are included in both of the strength 
function and the transition amplitude. Based on the 
analyses of the observed continuum and bound re- 
gion spectra, we make predictions to the future com- 
ing S~ hypernuclear production experiment at J- 
PARC. We find that we should observe E~ hypernu- 
clear bound state peak structures in (K~ , K + ) spec- 
tra on light nuclear targets such as 12 C and 27 A1 as 
far as the imaginary part of the S _ -nucleus optical 
potential is not large ( | We | < 3 MeV) and the exper- 
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imcntal resolution is good enough (AE < 2 McV). 

2. Green's function method and Local 
optimal Fermi averaging i-matrix 

The Green's function method in DWIA has been 
widely applied to analyse hypernuclear reactions. 
This method has an advantage that we can describe 
the continuum as well as bound state region on the 
same footing. In DWIA, the differential cross sec- 
tion reaction is obtained from the Fermi's golden 
rule [6], and the response function R(E) can be 
decmomposed into multipole components in the 
Green's function method [7], 

d2(T = Pk+E k+ 
dE K+ dQ K+ (2nh 2 ) 2 v K - { ' ' 1 ' 

R{E) = Y,\T f i\ 2 KE f -E i ), 
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where Vjc- is the incident K~ velocity, subscripts 
a and (3 stand for the quantum numbers of nucleoli 
and hyperon states, respectively, J is the total spin 
of hypernuclei, 6^ = 1 and for even and odd n, 
and </> Q (V) is the radial wave function of the tar- 
get nucleon. Dependence on the S~-nucleus optical 
potential Us appears through the Green's function 
G a f3a'i3'(E;r,r / ), which contains the hypernuclear 
Hamiltonian. The function j'jm is a radial part of 

( — )* (4-) 

the product of distorted waves Xk+ Xk~ evaluated 
in the eikonal approximation. We employ the tp ap- 
proximation for the imaginary part of distortion po- 
tential, Im[/if (r) = hvK^KN pij"), where <jkn is the 
isospin averaged cross sections with &nk- = 28.90 
mb and a N K + = 19.35 mb at P K - =1.65 GeV/c. For 
the real part, we adjust its strength to reproduce 
the total cross section data of K mesons [8]. The 
elementary i-matrix elements are usually assumed 
to be independent from the reaction point, and the 
Fermi averaging t-matrix squared are factorized [9] . 



For A and X productions, it is recently pointed out 
that on-shcll kinematics in the Fermi averaging (op- 
timal Fermi averaging, OFA) procedure roughly de- 
cide the shape of the QF spectrum [11] in the Green's 
function method with factrized t-matrix, and simi- 
lar procedure for t-matrix was proposed in Ref. [12]. 
In the Semi Classical Distorted Wave (SCDW) anal- 
yses [13], the local Fermi averaging of the elemen- 
tray cross section has been included. Here we would 
like to incorporate both of the above two ideas; we 
include the local optimal Fermi averaging t-matrix 
(LOFAt), t(r), in the integrand of the response func- 
tion Eq. (3). We define the LOFAt as, 

= Jdp N t( S ,t)p(p N )5HPf(r)-P?(r)) 

(r ' W ' 9j - Jdp N p(PN)S^(r)~Pt(r)) ' 

(6) 

where Pj*j (r) denote the four total momenta in the 
elementary initial and final two-body states. We 
adopt the Fermi distribution function for the target 
nucleon momentum distribution p(j>n) and param- 
eters are taken from [6,10]. In obtaining LOFAt, 
we define the i-th hadron single particle energy 
containing the nuclear and hypernuclear potential 
effects as, 

Ei(r) = y/pf + mf + 2miVi{r) ~ m t + ^-+V t (r) . 

(7) 

This treatment enables us to include the potential 
effects naturally through the effective mass m* 2 = 
mf + 2rriiVi{r), as adopted in transport models in 
high-energy heavy-ion collisions [14]. Consequently, 
the LOFAt has the dependence on the collision point 
r through hadron potentials Vi(r). 

3. Results 

In the calculation, we have assumed the one body 
Woods-Saxon type hyperon-nucleus optical poten- 
tial, U s (r) = (V S + iW?)f(r) + Vg (r), with /(r) = 
l/(l+exp((r-R)/d),R= r (^-l) 1/3 , d = 0.65 fm, 
ro = 1.1 fm, where V§{r) denotes 5 _ -core nucleus 
Coulomb potential. We assume the imaginary part 
of optical potential Wjf to be —1 McV, which sim- 
ulates the strength in the quark cluster model and 
the Nijmegen potential model D [15] estimations. 
We have adopted the elementary t-matrix, which is 
re-parameterized to fit the cross section and angle 
dependence for P K - < 3 GeV/c. 

Figure 1 shows the calculated results of S~ QF 
production spectra with potential depth of 14 MeV 
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Fig. 1. Calculated S — -hypernuclear production spectra in the QF region at p K - =1.65 GeV/c and 9 K + = 6 deg. on C, Al, Cu, 
Ag and Pb targets in comparison with data [5]. Solid lines show LOFAt + DWIA results with (V H , W^) = (— 14MeV, -lMeV), 
and dotted lines show the results without Kaon potential effects. In both of the calculations, the experimental resolution is 
assumed to be AE = 20 MeV (FWHM). 



in comparison with experimental data [5]. Calcu- 
lated curves reproduce the experimental data sys- 
tematically on heavy targets, Cu, Ag, and Pb, in 
the high p K + region, where the hypernuclear exci- 
tation is small. In the lower p K + region, other con- 
tributions have been known to be important [16,17], 
including heavy-meson production and its decay, 
K~N -> MY, M -> K+K- (M = </>,a ,/o) I 16 ! 
and the two-step strangeness exchange and produc- 
tion processes, K~N -> MY,MN -»• K+Y (M = 



it,r],p, ...) [17]. 

We underestimate the production spectra on 
lighter targets, 12 C and 27 A1. The underestimate of 
QF spectrum on 12 C target is a common feature in 
previous DWIA calculations [9,18]. In Ref. [13], it 
is discussed that this underestimate may be due to 
the center-of-mass effects: For electron scattering 
on a nucleus with mass number A, the center- 
of-mass correction in the shell model have been 
taken care of by a multiplicative factor F 1//2 = 



3 



cxp[q 2 /(AmiyAhuj)] for the form factor, where q is 
the momentum transfer. With fkj = 41A~ 1/3 McV, 
the factor [F 1/2 ] 2 amounts to 1.86, 1.43, 1.22, 1.16 
and 1.10 for 12 C, 27 Al, 63 Cu, 109 Ag and 208 Pb tar- 
gets, respectively, at q = 500 McV/c. In Ref. [19], 
we have adopted a different elementary i-matrix 
parameterization [17] and larger isospin-averaged 
KN cross sections, then we can roughly explain the 
target mass dependence, as shown with the dashed 
lines in Fig. 1. 
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Fig. 2. Calculated S production spectra on C target at 
P K - =1.65 (GeV/c). Calculated curves are shown for the po- 
tential depth of -24 MeV (dotted), -14 MeV (solid) and -4 
MeV (dashed). Thick lines show the results with multiplica- 
tive factors (1.7, 1.99 and 2.35 for V Q S = -24, -14, -4 MeV, 
respectively) introduced to fit the data, and thin lines show 
calculated results without these factors. 

Since our understanding is not complete and 
we have several ambiguities described above for 
the absolute yield in QF spectra, we introduce 
an adjustable multiplicative factor to fit the QF 
spectrum at low excitation energies (high K + mo- 
mentum region). In Fig. 2, we show the calculated 
S~ production spectra with 'Er potential depth of 
V E = -24, -14 and -4 MeV with multiplicative 
factors in comparison with data [5], The potential 
depth dependence is small in low p K + region, and 
attractive potential shifts the spectrum towards the 
high px+ direction slightly. When we adjust the 
multiplication factors as described above, calculated 
results reasonably well explain the QF spectrum. 
This means that we cannot determine the potential 
depth accurately from the QF spectrum shape. 

On the other hand, production spectrum is more 
sensitive to the potential depth at low excitation en- 
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Fig. 3. Potential depth dependence of the S — -hypernuclear 
production spectra in the bound state region at p,r=1.80 
GeV/c and 8 < 8 deg. on 12 C without(with) the experimental 
resolution (upper/lower panel). Dotted, solid and dashed 
lines show the results with E~-nucleus potential depths of 
24, 14 and 4 MeV, respectively. Experimental data are taken 
from Ref. [4], 

crgics, then there is a possibility that we can deter- 
mine the potential depth as discussed in Refs. [3]. 
In Fig. 3, we show the results of the potential depth 
dependence of the calculated S~ production spec- 
trum on 12 C target at p^- = 1-80 GeV/c multiplied 
by the factors described above in comparison with 
data [3]. 

Ikcda et al. evaluated the width of S~ hypernu- 
clear states ( 11 B+S _ ) to possible double A states 
based on the Nijmcgen model D potential as 1.2 



MeV and 0.5 MeV for [irQp 



3/2 



3 l/2j 



J=l 



and 



[7r0p 3 / 2 ® ^p] J_0 ' 2 , respectively [15]. In the upper 
panel of Fig. 3, we show the ideal production 
spectra on 12 C target without the energy resolu- 
tion folding. The imaginary part is assumed to be 
Wq = —1 MeV. These calculated spectra show 
that the S~ hypernuclear state widths are in good 
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Fig. 4. 



30 



25 



20 



2 15 



O 

a 

CO 
CO 
CO 

o 
« 



10 



1 ' W =-l'MeV 
W„=-3 MeV 



'C(K',K + ) 

P K =1.65 (GeV/c), 6 (deg.) 

V s =-14 (MeV) 
AE=2 (MeV) 




160 



140 



120 



■g 100 

O 80 

O 
W 

CO 
CO 
CO 

o 
M 

O 



60 



40 



20 



1 ' 1 ' 1 1 1 ' Wo = 
107 Ag(K- ; K + ) W ' = 


-lMeV ' 1 
-3 MeV 


_ P K =1.65 (GeV/c), 6 (deg.) i 


V z =-14(MeV) 
- AE=2(MeV) 


106 Pd + S- / 

1 / ■ 


f s 


a.,-7 o g9/2 -i" 






i .. J 1- t ' r i ' i 7 





-35 -30 -25 -20 -15 -10 -5 
-B^ (MeV) 



10 



■z, 
o 

o 
a 

CO 
CO 
CO 

O 
« 



> 

a 



o 
a 

CO 
CO 
CO 

O 
M 

O 



40 
35 
30 
25 
20 
15 
10 
5 



w = 

- 27 A1(K",K + ) Wo= 


-1 MeV — 

-3 MeV 


- P K =1.65 (GeV/c), 6 (deg 


.) - 


_ v q=-14 uvievj 
. AE=2(MeV) 


| ^Mg+S // 




d z // . 


Pe 




s s ..A-.. 


7'" " ' d 5/2~ 


. . . . j - -,- -i - ', r ', i 


1,1 



-25 



-20 



-15 



-10 



-5 
(MeV) 



10 



160 



140 



120 



100 



8 80 



60 



40 



20 



1 1 1 1 w l- 
208 pb(K - K+) Wo= 


-1 MeV 1 -7— 
-3 MeV 


P K =1.65 (GeV/c), 6 (deg 


M 


_ V S =-14 (MeV) 
_AE=2 (MeV) 

1 Se 


h,| V 

A j 207 ti + s- " 

m i ; 


f E / V 


*7 ■ / 


d z ,/V 


; ', ,' Ohn/2" 1 


_ .j _ - T - -~r ' , i 


1 



-40 



-30 



-20 -10 
-B^ (MeV) 



10 



production spectra at p K - = 1.65 GcV/c and K +=6 deg. on 12 C, 27 A1, 107 Ag and 208 Pb targets expected in 



the J-PARC experiment. We assume a Woods-Saxon potential with 14 MeV depth, and we show the results with imaginary 
parts of -1 MeV (solid) and -3 MeV (dotted). Experimental resolution is assumed to be AE = 2 MeV (FWHM). 



agreements with the estimates in Ref. [15]. In com- 
parison with experimental data, these spectra must 
be folded using a Gauss function. In the lower panel 
of Fig. 3, we show the results with an experimental 
resolution of AE = 12 MeV FWHM. Since the ex- 
perimental resolution is not enough to distinguish 
the bound state peaks and the statistics is low, we 
should compare the integrated yield in the bound 
state region. We find clear potential dependence in 
the bound state region, and with deep S~ potential 
(Us = —24 MeV) we may find a bump structure 
at around — £>- ~ — 12MeV even with this low res- 
olution. Comparison with the data suggests that 
the potential depth around 14 MeV is preferred 
in the present treatment, and if the statistics is 
high enough, it would be possible to determine the 



potential depth even with low resolution. 

The depth of the S~-nucleus potential has been 
already suggested to be around 14 MeV from the 
analysis of the (K~,K + ) spectrum in the bound 
state region [3]. In that analysis, the t-matrix ele- 
ment is evaluated under the frozen nucleon momen- 
tum approximation, where the kinematics is given 
with zero initial nucleoli momentum. In the present 
analysis, while the kinematics in the elementary pro- 
cess and the way to fix the absolute value are dif- 
ferent preferred potential depth is similar. This may 
be due to the fact that the excitation energy depen- 
dence of LOFAt is weak and smooth since the cov- 
ered K + momentum range is narrow in the bound 
state region. 

Now we find that S~-nucleus potential with 14 
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MeV depth well describes the spectra in the bound 
state and QF region for light nuclear targets, then 
it would be valuable to predict the peak structure 
which would be observed in the future coming J- 
PARC day-one experiment. Sensitivity of calculated 
spectra for the S~-nucleus potential is weak in the 
QF region because of the high momentum transfer 
q ~ 500 MeV/c, therefore it is difficult to extract 
precise potential information from the QF region, as 
shown in Fig. 2 

In Fig. 4, we show the calculated K + spectra 
in the bound state region of (K~,K + ) reactions 
on 12 C, 27 Al, 107 Ag and 208 Pb targets with a po- 
tential depth of Vq 1 = — 14 MeV, which explains 
the QF spectra and low resolution spectra in the 
bound region. We compare the results with Wjf = 
-1 MeV (solid lines) and = -3 MeV (dotted 
lines). We assume that the experimental resolution 
of AE = 2 MeV would be achieved. We find that 
bound state peaks are populated selectively due to 
high momentum transfer (q ~ 500 MeV/c) as in the 
A production spectra by (ir + , K + ) reactions (q ~ 
350 MeV/c), and these peaks can be identified in 
the high resolution experiment. 

In the Green's function method, target nucleon 
deep hole states are assumed to have large imagi- 
nary energies. Therefore, calculated results may be 
overestimating the spectra around the ground state 
due to the long Lorentzian tail from the deep hole 
states having finite contributions in this energy re- 
gion. This problem will be discussed in the future. 

4. Conclusion 

We have investigated the S~-nuclcus potential 
through cascade (S) hypernuclear production spec- 
tra by (K~ , K + ) reaction in the Green's function 
method [7] of the distorted wave impulse approxi- 
mation (DWIA) with the local optimal Fermi aver- 
aging i-matrix (LOFAt) [19]. The calculated spec- 
tra are in good agreement with the experimental 
data for heavy targets. With the multiplicative fac- 
tor adjusted to fit the spectra on light targets, we 
find that the calculated spectra well reproduce the 
observed spectrum [4] in shape and yield with S - - 
nuclcus potential depth around 14 MeV. This po- 
tential depth is consistent with those suggested in 
previous works [2,3,4]. 

While the dependence on the potential depth is 
small in the continuum region, it is clealy distin- 
guished in the bound region. Therefore, it would be 



possible to extract the S~-nucleus potential depth 
from the production yield in the bound state region 
when the statistics is high enough. Furthermore, the 
S - bound state peak structure can be found in the 
(K~,K + ) spectra on light target such as 12 C and 
27 A1, as far as the imaginary part is not very large 
(|W(f | < 3 MeV) and the experimental resolution 
is improved (AE ~ 2 MeV), as expected in the J- 
PARC experiment. We believe that our prediction 
would provide useful information in searching for 
the S~ nuclear bound states at J-PARC. 
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